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Abstract 

We present the theory for Multiple Energy X~ray Holography (MEXH) , using 
a multipole expansion for the scattered field. We find that light polarization 
plays a crucial role in the reconstruction of the image, and we suggest how 
to use it in order to eliminate aberration effects. The method we propose 
is alternative to the SWIFT method {scattered-wave-included Fourier trans- 
form), but has the advantage that no theoretical calculations are required to 

redefine the hologram. 
PACS 42.40.-i 
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I. INTRODUCTION. 



The idea of holography dates back to Gabor's original works0 and is essentially based on 
the interference of optical paths in crystal diffraction. The problem has been summarized 
very clearly by Sz5kei. Consider as a simple example, an emitter A at the origin and a 
single neighbouring atom A^^i at a position ri. In the process known as x-ray fluorescence 
holography (XFH)ifl, the atom A is excited and the emitted radiation either goes directly 
to the detector (photofilm) at a far distance R (the so called reference wave), or is firstly 
scattered by {object wave). The photofilm registers in this way the interference of the 
two optical paths (hologram), and gets in Fourier transform the desired information on the 
scatterer position ri. Equivalently, as proposed originally by Gaboifl, the last step can be 
accomplished by illuminating the photographic film with an incoming spherical waveffl (the 
decoding wave) , but numerical solutions that use a digital detector and suppress the decoding 
wave are nowadays preferable. 

A step forward in this direction has been given recently by the group of Fadley and 
Materliki, in using the emitter A as a detector. In this experiment (MEXH), the reference 
wave is furnished by a synchrotron radiation (SR) source, while the fluorescing atom A senses 
an electric field which is the sum of the direct wave and the scattered front (see Fig. |l]) due 
to its neighbouring atoms {object wave). 

Formally, MEXH can be visualized as a time reversal situation of XFH, in the sense that 
the detector is substituted by the source, which obviously does not modify the hologram. In 
spite of this, it presents many hidden advantages. 

First of all, the hologram x(k) is defined at different energies. In MEXH and in contrast 
to XFH, the photon energy is always above the fluorescence threshold of the emitter, thus 
allowing a better resolution of the image as ~ k~^. Secondly, the SR photon energy in 
MEXH can be freely varying, and this allows one to suppress unpleasant twin-image effects, 
once the integration is not limited to a sphere, but is done over a a 3-D volume in k spacei. 
A comparison between the two methods has been given recently by Len et a/.i and we refer 



to them for further details. 

In this report we develop the full vector theory for MEXH. The general theory is presented 
in section II. In Sec. Ill we make a multipole expansion of the previous formulae, and discuss 
their limit of applicability. Here is where aberration comes out, and we show in Sec. IV how 
to avoid it, playing with the polarization of the reference beam. Finally in Sec. V we 
confirm our previous findings presenting some numerical results on a real crystal structure. 
Particular attention is devoted to the disturbance induced by the scatterers at the large 
distance from the emitter. 



II. THE BASIC EQUATIONS. 

Let us fix our attention on the neighbouring atoms first. The external SR monochromatic 
field is represented by the normalized vector potential 

A(r, t) = ee^'^-'"e-*'^*+ c.c. (2.1) 

where k = uj/c and e is the polarization. This acts on each atom surrounding A as a 
disturbance Hj ~ — eA • p, causing a transition on the electronic states. The density and 
current involved in this process are 

p ~ -e ipnipm + {m ^ n) (2.2a) 
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(2.2b) 



where ipn is the initial core electron level, ipm is the final state, and the last term in (|2.2b| ) is 



the diamagnetic contribution. Each transition is then weighted by appropriate coefficients, 
which can be found by applying standard perturbation theory.1 Matching the density and 
the current density operators between states, the result is of the same form as in Eq. ( p.lD , 
i.e. 

p(r, t) = p(r; uj)e-''^* + c.c. (2.3a) 



and 



(2.3b) 



To find the two quantities on the r.h.s., we must pay attention that: 1) the atoms are 
centered at and not at the origin; 2) the wavefunctions ip depend on the relative position 
with respect to the nucleus 



r — r„ 



(2.4) 



and are otherwise independent of s. This means that the induced density (or current) for 
each atom, factorizes as a prefactor exp(ik ■ r^) times a term which is function of x in 
Eq. ( p.4|) and is formally independent of the atomic position. The total induced density (or 
current) is the sum over atoms. With this in mind the result is 



p{r;uj) = J2e^^-^^p{^) ; j(r; c.) = e^'^^'^ j(x) 



(2.5) 



where 



ip^^ * (x)V'^) (x) (m I e'-^'^pp \n) 
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In Eq. ( p.7| ) the dyadic over the momentum operator, p = —ih'V, means that it operates 
symmetrically both on the right and on the left, according to the square bracket in Eq. ( |2.2b[ ). 
In the first row in Eq. ( ^.T] ) one recognizes the diamagnetic contribution which does not 
depend on the virtual states m, while in the denominators we have defined 



- E^^) /h, 



and -E'^'*^ are energies of the excited and core electron level. It can be verified from ( p.6|) 
and ( |2.7| ) that Eq. (|2.3a|) and Eq. (|2.3b| ) satisfy the charge conservation V ■ j + p = 0. Once 



the charges and currents are given, the electric field at the emitter is calculated from the 
Maxwell equations. Outside the sources, this can be found from the vector potential only as 

E = ^VxB ; B = VxA (2.8) 

where 

A(r; u;) = - ff^iiv'; uj) dh' . 
c J \r — r I 

Changing the integration variable as in Eq. ( p.4|) , then using the definition in ( p.5|) , one gets 

E,^^{r;uj) = -J2e^-^4vx (2.9) 



/' jfc|r-rs-x'| \ 

^ xj(x') 

, r - r, - x' / 



The last equation is the object wave due to the scatterers. A similar result has been obtained 
recently by FondaS by computing the field through the vector and scalar potentials, but 
erroneously the latter was neglected from the start. Fonda's result, in the present notation, 
is Eobj = ikA, which is equivalent to dropping one term of the double vector product in 
Eq. ( [^.91 ). Further on, we show in Sec. IV how this correction becomes crucial for the 
polarization dependence of the reconstructed image, which is the main point of this paper. 
Adding from ( |2.1| ) the direct contribution of the reference wave, i.e.. 



E,ef(r; cj) = ifcee'''■^ (2.10) 



we get the total field 



Exot — Eref + Eobj (2-11) 

that acts as a perturbation on the fluorescing atom A. 



III. THE MULTIPOLE EXPANSION. 



The states of A involved in the transition are again core electron levels. This means that 
one needs Eq. ( 2.11 ) for r = x and |x| ^ d, with d of the order of the Bohr radius. In the 



same way, the integral in Eq. ( |2.9|) is confined to |x'| ^ d, i.e., in the volume where the 
current in ( |2.2b| ) is appreciably different from zero. Since is of the order of the lattice 
parameter, it follows that 



3> X, x' . 



This allows us to use for the propagator in Eq. ( p^) the so called plane wave approximation 
(PWA)0 

gi/c|r-rs— x'l gifc|r— r^l 
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e"*'^^-'' (3.1) 



with = —krg the scattered momentum. One gets 



Eobj(x; io) 
1 



^A;Vi'-''|[M,-(M,-f,)f, 
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[3(M,-f,)f , _ M,] 1 - ^ e*^-^^ (3.2) 



with 



M, = - / e-'''-"j(x')rfV. (3.3) 



Eq. ( p.2|) has been written for completeness. For kd <^ 1, exp(— iks-x') ^ 1 in Eq. (|3.3| ) and 
M, ^ p,§ with 

"xXxOrfV (3.4) 

the electric dipole moment, and p the density in Eq. (26). In this case the object field in 
( p.2[ ) is constant over the atom. In the two extreme limits kvg ^ 1 and krg <^ 1 one discovers 
in it the familiar result (see Ref. [ O]) of the radiative and the near static field in the dipole 



approximation. 



For x-rays the limit kvs ^ 1 is appropriate, but k ~ is required for good resolution. 
This means that under the sum in ( |3.2| ) only the first term is contributes, while in Eq. ( p.3|) 
one is led to attempt an expansion as 



The result is 



Eobj(x; u) = Y,k^e^^'"'\ [-f,x(f,xp)] 



+ [f,xm] - - [-f,x(f,xQJ] l^e**^-"-^ (3.5) 
where m is the magnetic dipole moment, 

(Qs)„ = EQa/3(k)^, (3.6a) 

/3=1 

and Qap is the electric quadrupole moment tensor 

Qap = J {Sx'^x'^ - r'^5^p) p(x') d^x'. (3.6b) 

To make the expansion consistent, we shall suppose in the following that the last two terms in 
Eq. ( |3.5| ) are small. Equally well we shall suppose that the same happens for the fluorescing 
atom A, such that the limit x = in Eq. (35) is appropriate. The yield is simply ~ |E(x = 
0)p and function of the direction k of the laser beam only. The hologram follows as 

X(k) = (lETotl' - lErefl') /|E,efr (3.7a) 
^ (E^f-Eobj + c.c.) (3.7b) 



and in the last line the self-hologram ~ E^j^j has been neglected. Using ( p.lO|) and ( p.4|) the 
result for Eq. ( |3.7b| ) is 

Akrs 



x{k) = -tkY^f,^e^^-^^ + c.c. (3.8) 
and where fg the scattering amplitude 



fs = (e*xfs) ■ (p X f J + e*xr,-m 

-^(e*xr,)-(Q,xf,). (3.9) 

Equation above shows that fs is complex and is polarization dependent. The quadrupole 
correction, causes a shifti, while the dipole terms give rise to a distortion of the image 
[aberration effects). In practical cases (see Ref. [ ^ for details) the shift is ~ 1/10 of the 
resolution and thus negligible. On the contrary the distortion is effective and is polarization 
dependent. 



IV. ANGULAR ANISOTROPIES. 

For a constant fs in Eq. (|3.8|) , the Fourier transform 

[/(r;fc) = i-|x(k)e-^'^"-dk, (4.1) 

with dk denoting the solid angle dQ\^, solves completely the problem of the atomic position. 
In this case@ 

^ ' ^ rs k\r-rs\ 

+^kY.ff-r^'^ML±Il\l (4.2) 
V Ts k\r + rs\ ^ ' 

and shows spherical illuminated spots with resolution ~ centered at the atoms and their 
twins. The twins do not represent a problem in MEXH. In fact, as suggested by Barton0, 
they can be mostly eliminated by integrating Eq. ( [4.1|) over the energy as 



Uir) = J U{r; k)e-''''' Pdk 



and we shall not discuss it. We shall concentrate on how to render fs constant. 

One way is to use the SWIFT method proposed by Saldin et a/.0. In it, the Fourier 
transform Eq. ( [4.1D is done not on the bare hologram x(k) as one gets from the data, but 
on the redefined quantity 



XswiFT(k; r) = ^ (4.3) 

Jr 



where 

fr = /s(f. f) 



and fs is defined in Eq. ( p.9|) . Using (|3.8|) one notes that as r in (^4. 1|) approaches r^, 



fs/fr ~ 1, which guarantees that both shifts and distortions are eliminated in Eq. (|4.1|) . 

The method applies for any polarization, but problems may arise when the ratio fs/ fr is 
~ 0/0. In addition the dipole moments and the quadrupole term are assumed to be known. 

Another way to render fs constant, would be instead to use the bare experimental quan- 
tity x(k), but to play with the polarization. 

Neglect for the moment the corrections in ( p.9|) and concentrate on the first term only. 
To make the life easier, here and in the following we shall suppose that p is in the direction 
of the field. From this we define 

~ |exf,|2 = sin^^^, (4.4) 

and fs = 1, if 

'"'^ (4.5) 



|f^xk| 

The choice made in Eq. ( [4.5|) implies the a priori knowledge of the position r^, but this can 
be roughly estimated through Eq. ( |4.1j ) and with an arbitrary polarization. The method we 
propose is the following. 

Let El, and €2 = kxei be two polarizations where the second is obtained by a rotation 
of 7r/2 over the direction k of the reference beam. Measure the yield for ei, €2 and for 

^3 = ^ (ei + £2) 

then call xi, X2 and xs the respective holograms. The third measurement collects interference 
effects and constructs the hologram X(q,/3) for any polarization 
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e = aei+pe2 ; |ap + |^|^ = 1 (4.6a) 



to be used in the experiment. For instance if both a and (3 are real one gets from (|3.7a|) the 
simple result 

XK/3) = (a - (3)iaxi - (3X2) + 2a/3x3- (4.6b) 
Taking for a and /3 the values 

Ir^xkl |rsXk| 

e = e^, thus = 1 as wanted. 

We show now with an example how, without corrections, distortion is effectively present. 
Suppose that is along the z axis and take 

z X k 



|z X k| 



in the (x, y) plane. Then leixr^p = 1, but 

|e2Xrsp = cos^ 'd 

where 

k = (sin-i^cosv?, sin sin cos-i?). 

Spherical spots are insured for ei, but not for €2 as we now show. 

Consider two paths: the first with r along and the second with r — in the (x, y) 
plane. For both paths the function V in Eq. ([4.1|) depends only on 



A simple calculation for the first path gives 



f^ii(C) = ^ 



1 - ^ j sin C + ^ COS C 



with a maximum ?7||(0) = 1/3, and vanishing at 
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For the second path U±{0) = f/||(0), but the shape is wider and vanishes at 

^ .757 

In this case then the illuminated spot is a highly anisotropic ellipse. The two axes of the 
ellipse are in the ratio 1.5/.757 ~ 2. 

V. NUMERICAL RESULTS AND CONCLUSIONS. 

A numerical calculation has been carried out for a Fe bcc structure, to see how the method 
works. With the emitter at the origin, we limit our attention on the nearest neighbour spot 
r* = where a = 2.87 Ais the lattice spacing for iron. The plane is the (110) one and 

the directions [110] and [001] as indicated in the figures. 

Fig. ^ is for polarization ei, here both orthogonal to f * and k as in Eq. ( [4.5|) . A spherical 
spot is clearly shown, in spite of the possible disturbances induced by other scatterers, and 
neglected in the analysis of the previous section. Fig. ^ refers to €2 polarization, gotten 
by rotating ei by 7r/2 around k as indicated previously. Here the intensity is lowered and 
the anisotropy is present indeed. Both figures are for reference beam energy E = 6.24 KeV. 
The beam energy cannot be arbitrary, but must be chosen to avoid out of phase overlap 
between an atom and its twin, a well known problem in single-energy holography!. In fact 
the sum in Eq. ( |3.8| ) is invariant for —>■ — r^, thus one gets (apart from a 2/c^ factor) the 
dimensionless quantity 

x(k)=E^^^T^cos(k.r,), (5.1) 

with /s < 1 for both polarizations given in Eq. ( ^4.4| ). The prefactor P = sin{krs)/krs selects 
the possible energies for the image reconstruction as sin(fcr*) ^ +1. Reversing the sign in 
defining the U function in Eq. ( |4.1| ), the values sm{kr*) ^ —1 are again possible, but we 
exclude them since are minima in the |ETotP intensity in Eq. ( p.7a| ). In the present case the 



most suitable values for the image reconstruction for r*-scatterer are E2 = 6.24 KeV and 
= 11.23 KeV (the first value = 1.25 KeV being too small). 

Figures and refer to this higher energy E^. Again one notices that the r*-atom is 



visualized only for ei polarization (Fig. [3a| ), while this spot vanishes using the polarization 

62 (Fig. 

To get closer to the experimental conditions, the x values in Eq. (|5.1| ) have been cal- 
culated on a 5° X 5°-grid over the whole solid angle, then averaged over an experimental 
window. A spline interpolation on these grid-data creates the averaged hologram to be used 
under the integral in Eq. ( [4.1|) . In all cases we find that the angular averaging is irrelevant 
as far as U{r) is concerned. The primary role is played by the integration over the modulus, 
i.e. over the energy. Here we use a 1-D gaussian convolution g ~ exp (-P/AA;^) on the 
grid data. The effect of using the gaussian is twofold, and well discussed in Ref. [ ||. On 
the one hand it is nothing but the energy spread of the primary beam. On the other, it 
acts as a low-pass filter to wash out all the distant scatterers, and at the same time avoids 
divergences in the reconstructed pattern. The first point is closely related to the dimension 
of the cluster, or mutatis mutandis, to the number of atoms of a bulk crystal we include under 
the sum in Eq. ( |3.8| ). We use a sphere of radius R = ax N and find that, for Ak ^ 7c/{aM) 
and N > M, the result is independent of the sphere's radius. 

There is however a divergence for large A^, which has more to do with our definition of 
U than with the physics itself. It can be avoided when the integration is not performed 
overthe whole solid angle as in Eq. ( [4.1|) , but only on a part of it. This is what is usually 
done by the people working in atomic holographyEi'0. Nevertheless we shall continue to use 
Eq. ( [4.1|) , and show how to get rid of this problem. 

Its origin can be figured out very simply by going to the scalar case, thus taking /s = 1 
constant. From Eq. (|4.2|) get the intensity at the origin 



which diverges for R —>■ 00. The same trend is of course present for |r| ^ 0, and this means 
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a very faint intensity in the whole pattern, except at the emitter. However, this divergence 
can be removed in the calculation by subtracting a C-constant, just as the reference wave 
is subtracted to get at the experimental x- The C value we use is fixed by taking constant 
U (0) as the sphere radius is increased from R = a to R = Na. This is a very good way to 
check how the result stabilizes with increasing number of scatterers, even if other choices 
are possible^!. 

We find much faster convergence at lower energy, as one would expect. The fourier 
transform (FT) U{r) is smoother and the illuminated spots are larger. 

Actually for E = E2 v/e do not even need to perform a gaussian convolution, but only 
a C -redefinition of the type discussed is required. Nevertheless an energy spread is used in 
all the figures we present. We fix a = ISkjk = 6.411 x 10~^, close to the Tegze and Feigeli 
value, this being the default value for R > 3a, and energy E^. We have checked, as just 
noted, that Figs. ^ and |2B| remain in practice the same using a larger a or setting Ak = 
(monocromatic beam). On the contrary, in Fig. ^ and Fig. ^ a gaussian convolution is 
essential. 

The numerical calculation is performed with a Fortran program. Once the averaged 
holograms Xi; X2 and xs are calculated on the 5° x 5° grid, the FT requires 2' and 30" of 



computer time for each reconstruction at E = 11.23 KeV. For Figs. ^ ^ the requirement 
reduces to ~ 1'. For the graphical part we use Mathematica. To get better contrast in the 
figures, we follow the standard procedure used in atomic holography^, namely plotting \U\'^ 
rather then \U\. Before doing this U has been redefined such that < U < 1 in the whole 
< X < V2a, < y < a (11, 0) plane, or in the relevant part of it. 

In conclusion, we can assert that the reconstruction of the image can be gotten with a rea- 
sonable effort from the experimental hologram. The energy spread of the nearly monochro- 
matic reference beam, prevents to see distant scatterers, but justifies to neglect the self- 
hologram contribution. An additional convolution of the data with a gaussian lowers the 
definition of the image, but makes the self-hologram even smaller. The polarization, together 
with the energy, plays an essential role for the observability of the atom. The polarization 
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' constraint ' remains effective even if, as in MEXH, the integration over energy is performed. 
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FIGURES 




FIG. 1. The picture of MEXH process. The interference pattern is due to the different optical 
paths of the direct beam (dashed) and the one scattered by the neighbouring atom Ni (soUd line). 
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a/^/J Vz'a 

[110] 

FIG. 2a. Theoretically generated atomic images for a Fe bcc crystal on the (110) plane. The 



X-ray energy is £^ = 6.24 KeV and the polarization ei is defined in the text. Other parameters are 



R = 3a (258 scatterers), C = .273, and Ak/k = .064 for the gaussian convolution. 
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FIG. 2b. As in Fig. 2a, but for the €2 polarization. 
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[110] 

FIG. 3a. Atomic images for E = 11.23 KeV and polarization ei. Here we have used i? = 4a 
(536 scatterers), C = .132, and the same Ak/k. 
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FIG. 3b. As in Fig. 3a, but for 62 polarization. 
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